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ABSTRACT 



We infer from different seismic observations the energy supplied per unit of time by turbulent convection to the acoustic 
modes of a Cen A (HD 128620), a star which is similar but not identical to the Sun. The inferred rates of energy supplied 
to the modes (i.e. mode excitation rates) are found to be significantly larger than in the Sun. They are compared with 
those computed with an excitation model that includes two sources of driving, the Reynolds stress contribution and 
the advection of entropy fluctuations. The model also uses a closure model, the Closure Model with Plumes (CMP 
hereafter), that takes the asymmetry between the up- and down-flows (i.e. the granules and plumes, respectively) 
into account. Different prescriptions for the eddy-time correlation function are also confronted to observational data. 
Calculations based on a Gaussian eddy-time correlation underestimate excitation rates compared with the values derived 
from observations for a Cen A. On the other hand, calculations based on a Lorentzian eddy-time correlation lie within 
the observational error bars. This confirms results obtained in the solar case. With respect to the helioseismic data, 
those obtained for a Cen A constitute an additional support for our model of excitation. We show that mode masses 
must be computed taking turbulent pressure into account. Finally, we emphasize the need for more accurate seismic 
measurements in order to discriminate, in the case of a Cen A, between the CMP closure model and the quasi-Normal 
Approximation as well as to confirm or not the need to include the excitation by the entropy fluctuations. 
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1. Introduction 

a Cen A is the most promising star after the Sun for con- 
straining the modelling of p-modc excitation by turbulent 
convection. Indeed, due to its proximity and its binarity, 
the fundamental parameters of a Cen A (effective tem- 
perature, luminosity, metallicity, gravity, radius) arc quite 
well known. For this reason this star and its companion 
(a Cen B) have been extensiv ely studied (see for instance 
the most recent modelling by iMiglio fc Montalbanl (|2005[ ) 
and the refer e nces t herein). As pointed out recently by 
I Samadi et al.l (|2007af ). the detection of p- modes and the 
measurement of their amplitudes as well as their mode line- 
widths (i.e. lifetime) from a Cen A enable to derive the 
rates at which energy is supplied to the acoustic modes for 
this star. These observational constraints can then be used 
to check models of p-mode excitation by turbulent convec- 
tion. 

Such comparisons have been first undertaken in the 
ca se of the Sun b y different authors (see the recent review 
by iHoudekl l2006f) . They enable to test different models of 
stochastic excitation of acoustic mode s as well as differen t 
models of turbulent convection (see eg. lSamadi et al.|[2006l ). 
Among those theoretica l prescriptions, we consider that of 
ISamadi fc Goupil (120011) w ith the improvements p roposed 
bv lSamadi et al.N2003l) andlBelkacem et all (|2006bD . It was 
shown bv lSamadi et alj that the way the eddy time 

correlation is modelled plays an important role on the effi- 
ciency of excitation. Indeed, calculations of the mode exci- 
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tation rates, V, that use a Lorentzian eddy-time correlation 
function better reproduce helio seismic data than those us- 
ing a Gaussian one. In addition, Belkacem et al.l (|2006bf ). in 
the case of the Sun, showed that excitation rates computed 
using an adapted closure model that takes the presence of 
plumes into account reproduce much better the solar obser- 
vations than the calcula tions based on the cla ssical Quasi- 
Normal Approximation (|Millionshchikov|| 1 94 If ) . 

An alternative theoretical model of the excitation 
of acoustic modes by turbulent convec tion proposed by 
IChaplin et alj (|2005f ) differs from that bv lSamadi fc Goupill 
( 20011 ) in several ways: it does not take the driving by 
the advection of the entropy fluctuations by the veloc- 
ity field into account. They only use the classical Quasi- 
Normal Approximation. More importantly, these authors 
claim that a Gaussian eddy-time correlation function re- 
produces better than a Lorentzian one the frequency de- 
pendence of mode excitation rates inferred from helioseis- 
mic data. However, they are led to introduce in their model 
a factor by which they multiply their formulation in order 
to reproduce the maximum of the solar mode excitation 
rates. 

a Cen A provides a second opportunity to test vari- 
ous assumptions in the modelling of the p-mode excita- 
tion: the amplitudes of the acoustic mod es detected in 
a Cen A were derived bv lButler et"aH (|2004[ ) using spectro- 
metric data. From those data, an estimate of the averaged 
mode line- widths has been first pro posed bvlBedding et al.l 
(|2004l ) and more recently updated in lKieldsen et al.l (|2005f ). 
Using a different met hod and data from the WIRE satellite, 
iFletcher et alj (|2006l ) proposed a new estimate of the aver- 



2 



Samadi R. et al. 



aged mode line- widths that diffe r significantly from the one 
derived bv lKjeldsen et alJ (|2005f ). Indeed, t he two data sets 
place the mode life ti me between 2.2 days (jKieldsen et al.l 
120051 ) and 3.9 days (|Fletcher et all I2006D . For compari- 
son, the averaged m ode life time derived for the Sun by 



Bedding et all 120041) in a similar way as for a Cen A by 
Kieldsen et alJ (120051) is about two days. 



ISamadi et alJ ( 2007a| ) have inferred from those sets of 
seismic constraints the p-modc excitation rates 7. They 
have found that they are significantly larger than those as- 
sociated with the solar p-modes. Furthermore, O 3 peaks in 
the frequency domain ~ 2.2 - 2.6 mHz while it peaks at 
the frequency v mal ~ 3.8 mHz in the case of the Sun. 

Although the spectroscopic characteristics (T e g = 
5810 K, log g = 4.305) of a Cen A are close to that of 
the Sun (T e g = 5780 K, log g = 4.438), the seismic signa- 
tures are quite different. Consequently, finding agreement 
between predicted and observed excitation rates would be 
a non-trivial result, providing additional support for the 
theory. 

A preliminary compar ison with theore t ical ca lculations 
obtained in the manner oflBelkacem et al.l (|2006a|) was car- 
ried out by ISamadi et alJ ( 2007al ). Discrepancies between 
the excitation rates inferred from the observations and the 
theoretical calculations were found. We update here this 
study by proceeding in a similar way as iRosenthal et alJ 
(|1999f ). Indeed, these authors have built a solar ID model 
where the surface layers arc taken directly from a fully 
compressible 3D hydrodynamical numerical model. We will 
refer here to such a ID model as a "patched" model. 
IRosenthal et alJ (| 19991 ) have obtained a much better agree- 
ment between observed and theoretical eigenfrequencies of 
the Sun computed for such a "patched" ID model than 
those obtained for a "standard" ID model based on the 
standard mixing-le ngth theory with no tur bulent pressure 
included. Following IRosenthal et~afl (|1999f ). we build here 
such a "patched" model to derive adiabatic mode radial 
eigen-displacements (£ r ) and mode inertia (I). We use them 
to compute the mode excitation rates, which we compare 
with excitation rates computed using £ r and / obtained 
with a "standard" ID model. 

The paper is organized as follows: in Sect. [5] we describe 
our procedure to compute the mode excitation rates for 
the specific case of a Ccn A. We then describe in Sect. [3] 
the way the mode excitation rates are inferred from seismic 
observations of a Cen A. In Sect.0] we compare theoretical 
calculations of 7 with those inferred from the seismic data 
obtained for a Cen A. Wc compare and explain in Sect[5]thc 
differences between a Cen A and the Sun. Finally, Sect. [5] 
and Sect. [7] are devoted to the discussion and conclusions, 
respectively. 

2. Modelling the excitation of p-modes 

2.1. General formulation 

The the oretical model of s t ochast ic excitati on is basically 
that of ISamadi fc Goupill (|200ll see also Samadi ct al. 

2005)) with the improvements of iBelkacem et al. 

200610), we thus recall here only some key features. The 



model takes two driving sources into account. The first 
one is related to the Reynolds stress tensor and as such 
represents a mechanical source of excitation. The second 
one is caused by the advection of the turbulent fluctuations 



of entropy by the turbulent motions (the so-called "entropy 
source ter m") and thus represents a thermal source o f 
excitation dGoldreich et al.lll99l IStein fc Nordlundll2001l) . 
The power fed into each radial mode, IP, is given by: 

y = jj{cl + cl) , (1) 

where C\ and C| are the turbulent Reynolds stress and 
entropy contributions, respectively and 

dm\tir\ 2 (2) 



is the mode inertia, £ r is the adiabatic radial mode dis- 
placement and M is the mass of the star. The expressions 
for C\ and G\ are given for a radial mode with frequency 
uj by 
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where f r = (d£ r /dr) 2 and g r is a function that involves 
the first and the second derivatives of £ r . E(k) is the spa- 
tial turbulent kinetic energy spectrum, E s (k) the spectrum 
associated with the entropy fluctuations, Xk the time cor- 
relation function of the eddies, a s = (dP/ds) p where s is 
the entropy, P the gas pressure, p the density, po the equi- 
librium density profile and loq the eigenfrequency. 



Finally, Si, 
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is the skew- 



ness and w the vertic al component of the velocity (sec 
IBelkacem et alJ l2006allbl . for details). Indeed, the expres- 
sion of Eq. © depends on the closure model used to ex- 
press the fourth-order moments involved in the theory in 
terms of the second order ones. The most commonly used 
closure model at the level of fourth-order moments is the 
quasi-normal approximation (QNA). Such an assumption 
leads to a vanishing skewness § w . However, in the solar 
case, the deviation from the QNA is due to the presence 
of turbulent plumes. By taking both the effect of the skew- 
ness introduced by the presence of two flo ws and the effect 
of turb ulence onto each flow into account. IBelkacem et al.l 
(l2006al ) thus proposed a new closure model that takes into 
account the presence of plumes, leading to a non- vanishing 
skewness, $ w , in Eq. ([3]). In the present work is then 
obtained directly from the 3D numerical model. 

Calculation of the mode e xcitation rates i s perform ed 
essentially in the manne r of Belkacem et al.l (|2006bf ) as 
detailed in ISamadi et al.l (|2007al ) in the specific case of 
a Cen A: all required quantities, except the mode eigen- 
functions £ r and mode inertia /, are directly obtained from 
a 3D simulation of the outer layers of a Cen A which char- 
acteristics are described in Sect. 12.21 below. 

The mode displacement £ r and mode inertia / must be 
computed from a global ID equilibrium model. We choose 
to study two such equilibrium models which arc described 
in Sect. 12.31 Finally, eigenfrequencies and eigenfunctions 
are computed using the adiabatic pulsation code ADIPLS 
([Christensen-Dalsgaard fc Berthomieulll99lf ). 
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2.2. The 3D hydrodynamical model of the outer layers of 
a Cen A 

We consider in this work the 3D hydrodyna mical model of 
the outer layers of a Cen A com puted by I Samadi et al.l 
(|2007al) using IStein fe Nordlundl l|l998l) 's code. The as- 
sumed micro-physics (e.g. the equation of s tate an d the 
opacity table) are detailed in I Samadi et al.l (|2007bf ). The 
hydrogen, helium and metal abundances are solar and the 
chemical mixture of the heavy elements is set according to 
iGrevesse fc Noel (fl~99l 's mixture. 

The 3D model associated with a Cen A has a hori- 
zontal size of 8.17 Mm x 8.17 Mm and a vertical size of 
4.31 Mm. Th e grid is 125 x 125 x 82. As pointed out by 
I Samadi et alj (|2007bf ). this spatial resolution is sufficient 
for the calculation of the p-mode excitation rates. The sim- 
ulation duration is 323 minutes while the acoustic depth of 
the simulation is 410 seconds and the characteristic eddy 
turn-over time is ~ 20 minutes (see Sect. [5]). The duration 
of the simulation then represents ~ 47 sound travels across 
the simulated domain and about 15 eddy turn-over times. 

The effective temperature T e g is adjusted to 5809 K ± 
15 in good agreement with the value T ft = 5810 K ± 50 
adopted bv lMiglio fc Montalbanl (|2005l ). The gravity is set 
to log g = 4.305 to match exactly the value (log g = 4.305 ± 
0.005) inferred from the precise measurements of the mass 
and the radius of the sta r (see the related references in 
iMielio fc Montalbanl[2005l) . 

2.3. ID models 
2.3.1. Standard model 

The first ID equilibrium model has the effective tempera- 
ture and the gravity of a Cen A and is built by imposing 
that for the temperature at the bottom of the 3D simulation 
box, the ID model has the same pressure and density as the 
3D simulation (see Fig. [TJ) . Hence, the 3D simulation is used 
to constrain this ID equilibrium model such that its interior 
structure is compatible with the second ID model described 
later on, in Sect. 12.3.21 (see also Fi g. []]). Convection in t he 
ID model is described according to iBohm-Vitensd (|1958[ )'s 
mixing-length local theory of convection (MLT) and tur- 
bulent pressure is ignored. Microscopic diffusion of helium 
and heavy el ements are treated accord ing to the simplified 
formalism of lMichaud fc Proffittl (119931) . We assume a solar 
abundance in order to be consistent with the 3D model. 

The mixing- length parameter, a, the age, the mass (M), 
the initial helium abundance (Yq) and the initial (Z/X)q 
ratio where X and Z are the hydrogen and metal mass frac- 
tions, respectively, are fitted such that the model simulta- 
neously reproduces the effective temperature of the star, its 
gravity, the solar composition and the temperature-pressure 
relation at the bottom of the 3D simulation box. The outer 
layers of this model, which matter here, then have the strat- 
ification given by a standard MLT model. This model will 
be referred to as standard hereafter. 

The matching results in a = 1.694. For comparison, 
the same matching performed for a solar 3D simulation 
results in a = 1.899. The mass of the standard model is 
M = 1.012 M© and the radius R = 1.1722 R Q . They arc 
sligthly smaller than that expected for this star, namely 
M = 1.105 ± 0.007 Mp and R = 1.224 ± 0.003 R Q 
(see iMiglio fc Montalbanl l^bOS) . This is because we have 



assumed a solar abundance for consistency with the 3D 
model. A global ID model with the iron to hydrogen abun- 
dance ([Fe/H]) of a Cen A (namely [Fe/H]=0.2), would 
have the expected mass and radius of the star. 

Slightly different R, M and [Fe/H] values might have 
some influence on the mode excitation rates (CP). In order 
to measure the effect of having an R, M and [Fc/H] different 
than required for a Cen A, we have computed two global 
models. The first model has an abundance [Fe/H]=0.2 and 
the second has a solar abundance. Both models have the 
effective temperature and gravity of a Cen A. In contrast 
with the "standard" model described above, we do not 
match these models with the 3D model. The model with 
[Fc/H]=0.2 has almost the radius and the mass expected 
for a Cen A while the second has almost the same R and 
M as the standard model investigated here. We find that 
y changes between the two models by less than ~ 5%; this 
is much smaller than the uncertainties associated with the 
observations. 

2.3.2. "Patched" model 

To take into account a more realistic description of 
the superadiabatic outer la yers, we have built, f ollowing 
iTrampedachl (|l997l see also ISamadi et all (|2007bD ). a ID 
global model in which the outer layers are replaced by the 
averaged 3D simulation (see Fig.[TJ). 

"Standard" and 3D models share the same micro- 
physics but mainly differ in the way convectivc motions and 
radiative transfer are treated. In the 3D model convectivc 
motions are treated by solving the Navier-Sokes equation 
while in the standard model convective motions are mod- 
elled according to the mixing-length model of convection 
and no turbulent pressure is included in the hydrostatic 
equation. Concerning the radiative transfer, in the standard 
model radiative transfer is gray and assumes the Eddington 
approximation. In the 3D model radiative transfer is ex- 
plicitly solved in LTE fo r four opacity bins (see details in 
IStein fc Nordlundl fl99l . 

The interior of the "patched" model is the same as in the 
standard model and does not include the turbulent pres- 
sure. At the bottom of the 3D simulation box, turbulent 
pressure is already negligible (~ 0.6 % of the total pressure). 
Then, neglecting it in the interior has negligible effects on 
the properties of the eigenfunctions considered here. This 
global model will be referred to as a patched model. Note 
that this patched model has the same total mass and a 
radius very close to that of the standard model, namely 
R = 1.1726 R . 

The stratifications in density of the patched and stan- 
dard models are compared in Fig. [TJ At the top of the 
convective region, we see that the density is lower in the 
patched model compared to the standard model. This is 
because the patched model includes turbulent pressure that 
provides additional support against gravity: accordingly, at 
a given total pressure (ftot)j the patched model has a lower 
gas pressure. Now, since T(P tot ) is the same with and with- 
out turbulent pressure, at a given temperautre the patched 
model has a lower density compared to the standard model. 

Because the treatment of photospheric radiative trans- 
fer is not (and in practice cannot be) identical between a 3D 
calculation and a ID model in the atmosphere, small dif- 
ferences in the stratifications of the very outer layers exist 
between the two models as seen in Fig. [TJ In any case these 
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differences do not play any significant role on the quantities 
which matter here, such as inertia. 

Note that some explanations about the differences seen 
between the outer layers o f 3D models and ID m odels 
have been proposed b y e.g. iNordlund fc SteirJ (|1999( ) and 
iRosenthal eFatl (|l999h . 




2.0x10* 1.5x10* 1.0x10* 5.0X10 3 

T [K] 



Fig. 1. Density as a function of the temperature. The solid 
line corresponds to the "patched" model and the dashed 
line to the "standard" model. The thick solid line is the 
part of the patched model obtained from the 3D simulation. 
The filled circle shows the position of the bottom of the 3D 
simulation box. 



3. Inferring the excitation rates from seismic 
constraints 

Mode excitation rates are derived from seismic observations 
according to the relation: 

?{v) = 2nMT{v/S ) 2 (7) 

where M = I/£ r (rh) is the mode mass evaluated at the 
layer = R+ h in the atmosphere where the mode is mea- 
sured in radial velocity, R the radius at the photosphere 
(i.e. at T = T e g), h the height above the photosphere, T 
the mode full width at half maximum (in v), v(rh, v) is the 
rms apparent velocity amplitude of the mode at the layer 
Th, v — ujq/2-k the mode frequency and Sq the visibility 
fac tor of the 1=0 mode. 

iKieldsen et alj (|2005f ) have derived the apparent am- 
plitude velocity spectrum, v(v) of the modes detected in 
a Cen A. However, their spectrum corresponds to ampli- 
tudes normalized to the mean of £=0 and 1 modes rather 
than to £=0. Furthermore, the y do not take into ac count 
the mode visibilities. Recently. iKieldsen et alj (|2008f ) have 
derived the (apparent) amplitudes of the modes, normalised 
to the £=0 modes and taking into account both the mode 
visibilities and limb-darkening effects. Finally, in order to 
derive the intrinsic mode amplitudes, we divide v(v) by 
Sq = 0.712, the visibility factor o f the 1=0 modes observed 
in velocity (|Kieldsen et alj|2008t ). 

For the mod e line-width, F, we us e the averaged val- 
ues pr ovided by IKieldsen et al.l (|2005[ ) and iFletcher et al.l 



Concerning mode masses, M, as discussed in Sect. [6j it 
is not trivial to determine the height h where the Doppler 
velocities are predominantly measured. As we do not pre- 
cisely know the height representative of the observations, 
we will next evaluate the mode masses - by default - at 
the optical depth r soonm — 0.013, which corresponds to 
the depth where the potassium (K) spectral line is formed 
(but see Sect. |5] for the discussion). This optical depth cor- 
responds to /i=470 km. 

Neither the standard nor the patched models do have 
exactly the radius and the mass expected for a Cen A (see 
Sect. 12.3. ip . However, this inconsistency has only negligi- 
ble effect on the mode mass M = Indeed, since the 
eigenmode displacement, £ r , is directly proportional to R, 
the mode inertia I scales as R 2 (see Eq. @). Accordingly, 
the ratio I/f£ is almost insensitive to a small change in R. 
Furthermore, we have checked that M is also insensitive to 
a small change in M. 



4. Comparison between observations and modelling 

We first compare theoretical calculations of 3 performed us- 
ing eigenfunctions computed with the patched equilibrium 
model with those computed using the standard equilibrium 
model (see Sect. 2.2). However, eigenmodes computed with 
those two models have not the same inertia and hence not 
the same mode masses M. Thus, we rather compare the ra- 
tios CP/M. As shown in Fig. [2j theoretical calculations that 
use the patched model lie well inside the observed domain 
of the seismic constraints. On the other hand, using the 
standard model leads to underestimated theoretical values 
compared to the two sets of seismic constraints. 

When comparing the integrands of the product TM - 
excitation power times mode mass - between the patched 
and standard equilibrium models, we find that they are 
quite similar. On the other hand, the mode masses M are 
quite different for the two equilibrium models in the domain 
1-3 mHz where the modes are mostly excited. This is due 
to the turbulent pressure which is present in the patched 
model and ignored in the standard model. At a given ra- 
dius in the super-adiabatic region, the patched model has 
a lower gas pressure and density (see Fig. [T]). As a conse- 
quence, inertia of the modes, which are confined within the 
super-adiabatic region, where the turbulent pressure has 
its maximum, are smaller for the patched model than for 
the standard model; accordingly, the ratios IP/M, which are 
inversely proportional to the squared mode mass M 2 , are 
about two times larger for the patched model. 

In Fig. [3l we compare two sets of calculations for a 
patched stellar model that assumes two different prescrip- 
tions for the eddy-time correlation function (xk) and two 
different closure models, namely the QNA and the CMP in 
the excitation model. The theoretical calculations based on 
a Lorentzian Xk and the CMP closure model lie inside the 
range allowed by the two sets of seismic constraints. The 
differences between calculations based on the CMP and on 
the QNA are found smaller than the differences between 
the two data sets. On the other hand, calculations based 
on a Gaussian Xk yields significantly underestimated val- 
ues compared to the seismic constraints. 

Note that I Samadi et alj (|2007af ) found a discrepancy be- 
tween theoretical calculations and observations. Part of this 
discrepancy was due to the fact that the horizontal size 
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of the simula tion box was set to tha t of the solar simu- 
lation used in lBelkacem et all (|2006bf K Indeed, the kinetic 
energy spectrum E involved in the expression for Sr and 
Ss in Eqs. (O and ([6]) must be normalised with respect to 
the horizontal size of the simulation box, as is done here. 
Further more, the mode inertia considered by I Samadi et all 
(|2007af) were computed for a standard MLT model (i.e. 
no turbulent pressure included) instead of using a patched 
model as is done here. As shown in Fig. [21 this results in 
an under-estimation of the ratio CP/M by a factor of about 
two (sec Fig. H}. 



Cen A 



a Cen A 




Fig. 2. Ratio of the rates CP at which energy is injected into 
p-modes to mode masses (M) for a Cen A. The dashed area 
represents the observed domain for CP/M = 2nY (v/So) 2 as 
a function of v. This domain is defined by merging the un- 
certainties associated with two independently derived val- 
ues of r and with the mode amplitudes v (Eq.[7J. The solid 
(resp. dashed) line corresponds to computed excitation 
rates with the eigenmodes obtained using the "patched" 
(resp. "standard") ID global model. All calculations here 
use the CMP and the Lorentzian function (LF) for the eddy 
time-correlation function \k in Eqs. [5] and [51 
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Fig. 3. Rates CP at which energy is injected into the p-modes 
of a Cen A. The dashed area has the same meaning as 
in Figure [21 The lines correspond to different theoretical 
calculations (all using a patched model) : the solid line uses 
the Lorentzian function (LF) and the CMP, the dashed line 
uses the LF and the QNA closure model, the dot-dashed 
line uses the Gaussian function (GF) and the CMP. 



If mode masses are evaluated at the photosphere (h = 0, 
T = Teg), we obtain 7 max = 15.9 ±8.0 x 10 15 [J/s]. In that 
case this is about ~ 4.4 ± 2 times larger than CP maXi0 . 

As seen in Fig. [H the frequency where 7 peaks is 
^2.4 mHz for a Cen A. For comparison, in the solar case, 
7 peaks around 3.8 mHz. Clearly, the modes in a Cen A 
are excited at lower frequency compared to the solar modes. 
Note also that the frequency domain where the derivation 
of CP is possible from the available seismic data is much 
smaller for a Cen A than for the Sun. This is obviously 
because the quality of the seismic data is much lower for 
a Cen A than for the Sun. 



5. Differences between a Cen A and the Sun 

5.1. Excitation rates 

Fig.0]compares the excitation rates, CP, inferred for a Cen A 
with those inferred from hclioscismic measurements ob- 
tained for the Sun. For a Cen A, excitation rates are 
obtained from the seismic measurements as explained in 
Sect. H 

For the Sun we con sider the helioscismic data studied 
bv iBaudin et~aTl (|2005h . Wc use here solar mode masses 
obtained with a patched model computed as for a Cen A 
in Sect. [21 Mode masses are evaluated for the optical depth 
t ~ 5 10 -4 since SOHO/GOLF obse rvations are b ased on 
the Na Dl and D2 spectral lines fsee lHoudekll2006f) . 

We find CP maXi0 ~ 3.5 ± 0.4 x 10 15 [J/s]. The excitation 
rates inferred for a Cen A with mode masses M evaluated 
at the optical depth associated with the potassium line (r ~ 
0.013) give CP max = 8.25 ± 1.0 x 10 15 [J/s]. This is about 
2.3 ± 0.3 times larger than CP max ,0- 
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Fig. 4. Excitation rates CP inferred from seismic data ac- 
cording to Eq. [71 Filled circles correspond to the helio- 
seismic constraints obtained bv IBaudin et al.l (|2005h . The 
dashed area represents the observed domain for the excita- 
tion rates derived for a Cen A. 
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5.2. Excitation rates as a function of depth 

Fig. [5] shows the integrand dV/dm of the excitation rates 
(Eq. I1I3I4P as a function of the temperature for the mode 
for which ft is maximum in the Sun and in a Cen A. 
The top panel shows the contribution of the Reynolds 
stress (d?-R./dm) and the bottom panel the contribution of 
the entropy fluctuations (df^/dm). Excitation due to the 
Reynolds stress is maximum where the rms value of veloc- 
ity, u, peaks. Excitation due to the entropy fluctuations is 
maximum where s, the rms value of the entropy fluctua- 
tions, peaks. Fig. [5] shows that the excitation rate is larger 
for a Cen A than for the Sun and occurs over a slightly 
more extended region in a Cen A than in the Sun. 

The excitation due to the entropy fluctuations occurs in 
a more shallow region compared to the Reynolds stress. For 
a Cen A, the relative contribution of the entropy fluctua- 
tions to the total excitation is ~ 18%, which is similar to 
the Sun (~ 15 %). Hence, in both cases the excitation due 
to the Reynolds stress remains the dominant contribution. 

5.3. Differences in the characteristic properties of convection 

To summarize, we find that T is significantly larger in 
a Cen A than in the Sun. Furthermore, ft peaks at lower 
frequency. As shown below, all these seismic differences can 
be attributed to differences in the characteristic properties 
of convection between a Cen A and the Sun. 

5.3.1. Why 7 is larger for a Cen A ? 

At a given layer, the power supplied to the modes by the 
Reynolds stress is proportional - per unit mass - to 
po u 3 A 4 where A is the char acteristic eddy size and u the 
rms value of the velocity (see lSamadi fc Goupij[200lh . The 
flux of the kinetic energy, Fkin, is proportional to pou 3 . 
Hence, the larger F^[ n or A, the larger the driving by the 
Reynolds stress. 

The power supplied to the modes by the so-called en- 
tropy source term is proportional - per unit mass 
to po u* 3 A 4 ft 2 /(ta luo) 2 where w is the mode frequency, 
t\ rsj K/u is the characteristic eddy turn over time, finally 
ft oc -Fconv/-Fkin where F conv oc w a s s is the c onvective flux 
and s is the rms of entropy fluctuations (see ISamadi et alj 
|2006| ). Note that the higher ft, the higher the relative con- 
tribution of the entropy source to the excitation. Note also 
that the driving is max imum for mode frequency (see, e.g., 
ISamadi fc Goupilll2001h 

c^o ~ 2tt/t a . (8) 

Hence, at the mode frequency ujq ~ 27t/ta, the larger the 
ratio fconvZ-Fkin the larger the relative contribution of the 
entropy source term to the total excitation rate. 

As a summary, for both Reynolds stress and entropy 
contributions, the larger the characteristic scale length (A) 
or the larger the kinetic energy (Fkin), the larger the ex- 
citation. Furthermore, the larger ft, the larger the relative 
contribution of the entropy source term to the excitation. 
We study below the differences in A, Fki n and ft between 
the Sun and a Cen A. 

Kinetic energy flux (F^): 
The maximum in u is up to ~ 10 % larger in the 3D sim- 
ulation associated with a Cen A than in the solar one. 
However, the differences in the flux of kinetic energy, F^ n , 
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Fig. 5. Top: The integrand dVji/dm (Eq. [T]) associated 
with the contribution of the Reynolds stress to the excita- 
tion is plotted as a function of the horizontally and tem- 
porally averaged temperature in the simulation box for the 
mode for which ft is maximum. The solid line corresponds 
to the 3D simulation associated with a Cen A and the 
dashed line to the one associated with the Sun. Bottom: 
as the top panel for dfts/dm, the integrand associated with 
the contribution due to the entropy fluctuations. 



between the 3D simulations associated with a Cen A and 
the solar one are small (< 10 %). This small effect on Fkin 
despite its cubic dependence on u is due the lower po for a 
layer with the same average T in the simulation for a Cen A 
as compared to the simulation for the Sun. The lower po in 
turn is a consequence of the lower surface gravity of a Cen A 
compared to the Sun. 

Relative contribution of the entropy source term: 
We also find that s is ~ 25 % larger in the 3D simu- 
lation associated with a Cen A. However, the convective 
flux, F conv oc w a s s, in the 3D simulation associated with 
a Cen A is very close to that of the solar simulation. This 
is not surprising since the two stars have almost the same 
effective temperature. Furthermore, as pointed out above, 
the differences in i^km between a Cen A and the Sun are 
small. As a consequence ft oc F conv /F^ n does not differ 
between a Cen A and the Sun. This explains why the con- 
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tribution of the entropy term relative to the Reynolds stress 
is similar between a Cen A and the Sun. 

Characteristic scale length (A) : 
Fig. [6] shows the kinetic energy spectrum E as a func- 
tion of the horizontal wavenumber k and the scale length 
Afc = 2n/k for the layer where u is maximum. As seen 
in Fig. [5J for the 3D simulation associated with a Cen A, 
E is maximum at a larger scale length compared to the 
solar simulation. Then, the eddies have a larger character- 
istic scale length in a Cen A than in the Sun. This explains 
why the excitation of p-modes is significantly stronger for 
a Cen A than for the Sun. 

Note that since the number of grid points is the same 
for both simulations, the a Cen A simulation has a larger 
physical grid size, thus a smaller maximum wavenumber, 
hence the cut-off in the spectrum occurs at a lower k. This 
explains the earlier drop-off of E(k) for a Cen A in Fig. [6] 
For that reason the high wavenumber part (beyond a k 
value of about 15 Mm -1 , or a A& less than 0.4 Mm) should 
not be compared directly. On the other hand, the scaling 
chosen in Fig. [7J allows a direct comparison. 




Fig. 7. The kinetic energy spectrum, E, as a function of 
k Hp (lower axis) and A^/Hp (upper axis) for the layer 
where u is maximum. The lines used have the same meaning 
as in Fig. [6l 
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Fig. 6. The kinetic energy spectrum, E, as a function of the 
horizontal wavenumber k (lower axis) and the scale length 
Afc = 2-k /k (upper axis) for the layer where u is maximum. 
The solid line corresponds to the 3D simulation associated 
with a Cen A and the dashed line to the one associated 
with the Sun. 



We point out that the characteristic scale length, A, 
scales as the pressure scale height. Indeed, we have plotted 
in Fig. [7J the kinetic energy spectrum, E, as a function of 
kHp where H p is the pressure height at the layer where u 
is maximum. Except at small scale lengths, we see that the 
fc-dependency of the spectrum is almost the same between 
the simulation associated with a Cen A and the solar one. 

The ratio between H* ccnA and Hf ( H« ccnA /Hf ~ 
1.38) is very close to the ratio go/ g a ccnA ( — 1.36). This 
is obviously related to the fact that H p = P/pg cx T/g. 
Accordingly, since f a A 4 (sec above), we then have 
3>/3> cx (A/A ) 4 (x cenA /i?"®) 4 cx ( 5q cc„a/5o) 4 ~ 
3.4. For comparison excitation rates computed for a Cen A 
are two times larger than in the Sun. 



5.3.2. Why CP peaks at lower frequency ? 

The characteristic eddy turn-over time can be estimated as 
the quantity f ~ L^/u where Lh is the horizontal extend 
of the 3D model and u the velocity at a given layer. At the 
layer where u is maximum, we find that f, evaluated at the 
layer where u peaks, is larger in a Cen A (~ 30 minutes) 
than in the Sun (~ 23 minutes). This explains the fact that 
for a Cen A IP peaks at lower frequency than in the Sun 
(lo ~ 27r/f, cf. Eq.E). 

Note that both u and A are larger for a Cen A than for 
the Sun. However, the net result is a larger f for a Cen A. 

5.3.3. Interpretation 

The differences in characteristics of convection between 
a Cen A and the Sun can be understood as follows: as 
seen in Sect 15. 3.11 the characteristic size A is mainly con- 
trolled by Hp on T/g (for a given composition). The surface 
gravity for a Cen A is ~ 35 % times smaller compared to 
the Sun while the effective temperature is very similar to 
that of the Sun. Consequently, A is larger than in the Sun. 
Furthermore because of the lower gravity, the density at 
the photosphere is smaller than in the Sun. Hence, in or- 
der to transport by convection the same amount of energy 
per unit surface area, the convective cells must have larger 
speed (it). 

6. Discussion 

6.1. EfFect of chemical composition 

The star a Cen A has an iron to hydroge n abundance of 
[Fe/H]=0.2 fsee lMiglio fc Montalbanll2005D . The 3D simu- 
lation considered here has a solar abundance. Preliminary 
work tends to show that, at the given effective temperature, 
a 3D simulation with a metal abundance 10 times smaller 
than the solar one results in mode excitation rates ~ 2 
times smaller. This can be understood as follows: the ra- 
diative flux is larger for a low metallicity than in a medium 
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with a solar mctallicity. In that case, to transport the same 
amount of energy, convection is less vigourous (i.e. lower 
flux of kinetic energy, Fkin), leading to a lower efficiency of 
the driving. If we extrapolate this preliminary result, we can 
expect that mode excitation rates ought to increase with 
[Fc/H]. A quantitative estimate of the expected increase 
must however be performed, in particular for a Cen A, 
which will require to compute a 3D simulation with a non- 
solar abundance representative of the surface layers of the 
star (in progress). 

6.2. Estimation of mode mass 

Mode masses must be evaluated at the layer in the atmo- 
sphere where the acoustic modes are predominantly mea- 
sured. The result of a comparison with seismic constraints 
significantly depends on the effective heights h where mode 
masses are evaluated. Indeed, we plot in Fig. [8] the ratio 
T/M for mode masses evaluated at different heights h in 
the atmosphere, namely from h = (the photosphere) up 
to the top of the simulated domain (h ~ 1000 km). This ra- 
tio is compared to the quantity 2-kY (v/So) 2 obtained from 
the seismic constraints (Eq. [7]). For h > 600 km (i.e. for 
optical depth < 0.005), the ratio T/M. is outside the obser- 
vational domain. 



lar seismic measurements it is then more difficult than for 
helioseismic observations to es timate the effective height 
h (for the sol ar case see, e.g . . iBaudin et al.ll2005l ). A re- 
cent work by IKieldsen et all (|2008f ) allows us to estimate 
the value for an effective h. Indeed, the authors have found 
that solar modes measured with the UCLES spectrograph 
have amplitudes slightly smaller than those measured by 
the BiSON network. The instruments of the BiSON net- 
work use the potassium (K) spectral line which is forme d 
at an optical depth t 5 qq nm ~ 0.013 fsee iHoudekl I2006T ). 
IKieldsen et al.l (|2008f) 's results then suggest that acoustic 
modes measured by UCLES arc measured at an effective 
height (h) slightly below the formation depth of the K 
line, i.e. at optical depth slightly above T500 nm — 0.013. 
Accordingly, we have here evaluated the mode masses at 
that optical depth, which corresponds to h~A7Q km. 

A more rigorous approach would be to compute an ef- 
fective mode mass by weighting appropriately the different 
mode masses associated with the different spectral lines 
that contribute to the seismic measure. In order to infer 
accurate mode excitation rates from the seismic data of 
a Cen A, the mode masses representative for the observa- 
tion technique and the spectral lines of a Cen A must be 
derived. This is out of scope here. 



Cen A 



1 00 1 




Fig. 8. Ratio of the rates 7 at which energy is injected into 
the p-modes to the mode masses (M) for a Cen A. The 
dashed area represents the observed domain for !P/M = 
2ttT(v/So) 2 as a function of v (see Sect. [3]). The solid lines 
correspond to the ratio V/M where the excitation rates, 
are calculated according to Eq. Q]and the mode masses, 
M, are evaluated at different heights h above the photo- 
sphere. The lower curve corresponds to the photosphere 
(h = 0) and the upper curve to the top of the atmosphere 
(h = 1000 km). The step in h is 200 km. 



Seismic observations of a Cen A were performed us- 
ing UCLES and UVES spectrographs. UVES and UCLES 
use a similar technique to measure the acoustic modes (T. 
Bedding, private communication). As other spectrographs 
dedicated to stellar seismic measurements, the UCLES in- 
strument uses a large number of spectral lines in order to 
reach a high enough signal to noise ratio. In the case of stel- 



7. Conclusions 

Theoretical estimations for the energy supplied per unit 
of time by turbulent convection (O 3 ) to a Cen A acoustic 
modes have been compared to values obtained from obser- 
vations. This allows us to draw the following conclusions: 



7.1. Differences with the Sun 

Although a Cen A has an effective temperature very close 
to that of the Sun, we find here that the p-mode excita- 
tion rates CP inferred from the seismic constraints obtained 
for a Cen A are about two times larger than in the Sun. 
These differences are attributed to the fact that the eddies 
in a Cen A have a larger characteristic size (A) than in the 
Sun. This is related to the fact that a Cen A has a smaller 
surface gravity. 

Furthermore, the p-mode excitation rates for a Cen A 
are maximum at lower frequencies than in the Sun. This 
behaviour is related to the fact that the eddies have a longer 
turn-over time as a result of a larger A. 

The seismic characteristics of the p-modes detected in 
a Cen A significantly differ from that of the Sun. They 
can therefore provide additional constraints on the model 
of stochastic excitation. 



7.2. Inferred versus modelled excitation rates 

Our modelling gives rise to excitation rates within the er- 
ror bars associated with the observational constraints. We 
stress that this modelling was undertaken for a Cen A in- 
dependently from the solar case, i.e. without using a for- 
mulation fitted on the helioseismic data as it is the case , 
for instance, in the case of the Sun in Chaplin et al.l (|2005h 
or in the case of a Cen A in IHoudekl (I2003 T The seismic 
constraints from a Cen A then provide a clear validation of 
the basic underlying physical assumptions included in the 
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theoretical model of stochastic excitation, at least for stars 
not too different from the Sun. 



7.3. Constraints on the description of turbulence: eddy-time 
correlation 

We find that our theoretical estimations of 3 , which as- 
sume a Lorcntzian eddy-time correlation function (xk) 
and the Closur e Model with Plumes (CMP) proposed by 
iBelkacem et al.l (|2006al ). lie in the observed domain. On 
the other hand, when a Gaussian function is chosen for Xk, 
7 is significantly underestimated. The comparison with the 
seismic data for a Cen A confi r ms th e results for the solar 
case obtained by I Samadi et al.l (|2003l) that Xk s i gnific antly 
departs from a Gaussian. As in ISamadi et al.l (|2003D . we 
attribute the departure of Xk from a Gaussian to diving 
plumes (i.e. down-flows), which are more turbulent than 
granules (i.e. the up-flows). This result confirms that a 
Lorentzian function is a more adequate description for the 
eddy-time correlation than a Gaussian. 



7.4. Constraints on the modelling of turbulent convection in 
the equilibrium stellar model 

Calculations involving cigenfunctions computed on the ba- 
sis of a global ID model that includes a realistic descrip- 
tion of the outer layers of the star (taken from 3D sim- 
ulations) reproduce much better (see Fig. [2]) the seismic 
data than calculations that use cigenfunctions computed 
with a standard stellar model built with the mixing-length 
theory (MLT) and ignoring turbulent pressure. This is be- 
cause a model that includes turbulent pressure results in 
larger mode masses M than a model which ignores turbu- 
lent pressure. This can be understood as follows. Within 
the supcr-adiabatic region, a model that includes turbulent 
pressure provides an additional support against gravity and 
hence has a lower gas pressure and density (see Fig.[T]) than 
a model that does not include turbulent pressure. As a con- 
sequence, mode inertia (and hence mode masses) are then 
larger in a model that includes turbulent pressure. 

These conclusion are similar to that obt ained in the Sun. 
Indeed, the mode masses considered by IBelkacem et al.l 
12006b!) in the case of the Sun were obtained with a ID 
model computed using iGouehl (|1977l) 's non-local mixing- 
length formulation of convection. The model thus includes 
turbulent pressure. We do not observe significant differences 
between excitation rates obtained with this non-local model 
and those obtained with a "patched" solar computed as de- 
scribed here in the case of a Cen A. On the other hand, 
excitation rates computed with mode masses obtained with 
a "standard" solar model (that is with no turbulent pres- 
sure included) or with a model in which turbulent pressure 
is included according to the mixing-length theory under- 
estimate significantly the helioseismic constraints. 

These results tell that one must compute mode masses 
from ID models that include turbulent pressure using a 3D 
hydrodynamical model or using a non-local description of 
convection. 



7.5. Need for improved data sets 



As shown by ISamadi et al.l (|2003h in the case of the Sun, 
contribution of the entropy fluctuations to the excitation 



cannot be neglected. Furthermore, recently. IBelkacem et al.l 
(|2006bf ) have shown that theoretical calculations based on 
the CMP result in a better agreement with the helioseis- 
mic constraints than those based on the Quasi-Normal 
Approximation (QNA). 

However, in the case of a Cen A, differences between 
theoretical calculations that use the CMP and those based 
on the QNA (see Fig. [3]) as well as differences between calcu- 
lations including driving by entropy fluctuations and those 
that do not include it (not shown) , are of the same order as 
the observational uncertainties associated with the two data 
sets. The present seismic constraints therefore arc unable to 
discriminate between these assumptions. This emphasizes 
the need for more accurate seismic data for a Cen A. 
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